Vortex dynamics in NbTi films at high frequency and high DC magnetic fields

We report on the characterization of NbTi films at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼ 11 GHz and in DC magnetic fields up to 4 T, performed by means of the coplanar waveguide resonator technique, providing quantitative information about the penetration depth, the complex impedance, and the vortex-motion-induced complex resistivity. This kind of characterization is essential for the development of radiofrequency cavity technology. To access the vortex-pinning parameters, the complex impedance was analyzed within the formalism of the Campbell penetration depth. Measurements in this frequency range allowed us to determine the complete set of vortex-pinning parameters and the flux flow resistivity, both analyzed and discussed in the framework of high-frequency vortex dynamics models. The analysis also benefits from the comparison with results obtained by a dielectric-loaded resonator technique on similar samples and by other ancillary structural and electromagnetic characterization techniques that provide us with a comprehensive picture of the material. It turns out that the normalized flux flow resistivity follows remarkably well the trend predicted by the time dependent Ginzburg-Landau theory, while the pinning constant exhibits a decreasing trend with the field which points to a collective pinning regime.


Results and discussion
CPWRs were produced by the process described in the "Methods" section. Two devices are analyzed in this work, representative of different deposition processes. In the two devices, hereafter labelled as #1 and #2, the central strip of the CPWR has the following dimensions: length-width-thickness = 8.6 mm-282 μm-1.43 μm (#1), and 8.5 mm-244 μm-2.4 μm (#2). With these geometries, the resonance frequency of both the CPWRs sets around 11 GHz. A sketch of the CPWR (not to scale), showing the directions of the rf magnetic field and of the applied DC field, is shown in Fig. 1.

London penetration depth.
Resonance curves of the CPWRs were acquired and fitted by the procedure described in the "Methods" section. Examples of curves in zero DC field and at different temperatures are shown in Fig. 1. The inset shows the temperature dependence of the fitting parameters-the resonance frequency f 0 and the inverse of the unloaded quality factor Q 0 .
To obtain an estimation of the absolute value of the London penetration depth, L , we investigated the resonance frequency dependence on temperature, in the absence of DC magnetic field. The resonance frequency depends on L (T) through the inductance of the line, and on the permittivity of the substrate ε r (T) through the capacitance of the line. In fact, in the theory of distributed element transmission lines for an half-wavelength resonator, f 0 is given by www.nature.com/scientificreports/ where l is the length of the resonator, L g l is the geometrical inductance per unit length, L k l is the kinetic inductance per unit length, and C l is the capacitance per unit length. The way to explicitly express C l , L g l , and L k l as a function of L (T) , of the geometrical parameters of the CPWR and of the permittivity of the substrate ε r (T) is rather intricate and fully described in Refs. [7][8][9] . As for the temperature dependencies of L (T) and ε r (T) , we adopted parametric functions allowing easy data fitting. Going into details, for ε r (T) of quartz, literature shows that at microwave frequencies it faintly depends on temperature, within the analyzed range 10 , and that its absolute values range from 4.4 to 4.6, depending on the orientation. For the sake of generality, in the fit we assume a second-order polynomial temperature dependence for ε r (T) , even though we expect small linear and quadratic coefficients. Then, the way it enters in Eq. (1) depends on the geometry of the CPWR 11 .
The expression to be used for L (T) needs a deeper insight. With the aim to avoid forcing a priori a given trend, we used a very general function, that well adapts to different scenarios: From a phenomenological point of view, with γ = 4 Eq. (2) corresponds to the two fluids result, while γ ≈ 2 was typically obtained for d-wave superconductors (e.g. YBa 2 Cu 3 O 7−x 12 ). On a more fundamental basis, Eq. (2) also adapts to BCS-based calculations: the weak coupling clean limit is well approximated by γ = 3 − T/T c 13 (as successfully used e.g. with MgB 2 11 ). However, Nb, Nb 3 Sn and other Nb-based alloys are better described within the strong coupling regime. Curiously, this is well approximated by the phenomenological two fluids formula, with γ slightly lower than 4 for the clean case and γ slightly larger than 4 for the dirty case 14 . Therefore, in data fitting we adopted Eq. (2), with γ as a fitting parameter.
For the practical fitting procedure, we operated with the normalized resonance frequency that, through Eq. (1), can be written as: where T 0 is the lowest measured temperature. Taking into account all the assumptions discussed above, we fit the f 0 (T)/f 0 (T 0 ) data with L (0), T c , γ , and the second-order coefficients for ε r (T) as fitting parameters. Further details of this fitting procedure are reported in Ref. 9 , where the way to choose a given temperature range for the fit and then to extend the determination of L to the whole range, and the way to determine the surface resistance R s are described. An example is given in Fig. 2a for CPWR#1. Results are shown in Fig. 2b,c, in terms of L (T) and R s (T).
The R s (T) curves show a very good agreement between the two CPWRs. As for the penetration depth, it should be noted that, due to the thickness d of the films, the condition passes from the bulk limit at low temperatures, where 2 L < d , to the thin film limit at higher temperatures, where 2 L > d . In the latter case, www.nature.com/scientificreports/ the effective penetration length is the Pearl length, 2 2 L /d . To take into account both these limits, we used the phenomenological interpolation expression = L coth(2 2 L /d) to analyze the response of each CPWR, with its specific thickness. Then, we extracted the London penetration depth, which is characteristic of the material and independent of the CPWR dimensions, that is the one shown in Fig. 2b: the two curves reported, corresponding to the two CPWRs, define the L uncertainty for what concerns sample-to-sample repeatability. The resulting London penetration depth for T → 0 , as obtained on the basis of Eq. (2), is L (0) = (290 ± 40) nm for CPWR #1 and L (0) = (300 ± 30) nm for CPWR #2, and γ ≈ 4.1 for both, i.e. consistent with the strong coupling dirty case. These estimations of L (0) are in accordance with results reported in Ref. 15 (thin films, experimental) and Ref. 16 (theoretical calculations) for NbTi with similar Ti contents.
Effects of a DC magnetic field. The measurements as a function of temperature were repeated for different DC magnetic fields applied perpendicular to the CPWR surface. From these measurements, the values of T c (H dc ) can be extracted (as the onset of a superconductive response over the noise), and the upper critical field B c2 (T) can be evaluated as the applied µ 0 H dc field at which T = T c (H dc ) . Figure 3 shows the data measured with both the CPWRs (solid symbols), as a function of the reduced temperature, i.e. the temperature normalized by the zero-field critical temperature T c0 . CPWR data are consistent with B c2 (T) data obtained from DC resistance measurements R(T) of the Hall bars (red crosses), by the criterion of selecting the temperature at which the resistance reaches 90% of the normal state value, R n . Note that all the curves, from the two CPWRs and from DC resistance, show a downward B c2 (T) curvature, different from what expected by the standard Ginzburg-Landau (GL) theory. However, the curvature can still be understood within the GL theory if the granularity of the films is considered. In fact, approaching T c , the coherence length ξ GL divergence is cut off at a typical length scale corresponding to the grain size and/or void spacing. Above this point, the B c2 curve should assume a (1 − T/T c ) 1/2 behavior 17 . Within this frame, it was also proposed that a B c2 (T) dependence not as sharp as a square root variation can be justified by a distribution of grain sizes 17 . The presence of a granular structure supporting this hypothesis (even if with well-connected grains) is visible in our films by FESEM and can be inferred from the DC resistance normal-to-superconductor transition, with moderately large width of 0.4-0.5 K.
In Fig. 3 we show the fit (solid lines) of the CPWRs data with a generic phenomenological function B c2,0 [1 − (T/T c ) β ] (values of β of 4.1 and 6.6), extrapolated to lower temperatures (dotted lines) to check the matching with two further points (open symbols), obtained by the TDGL analysis of the ρ ff (T) data, as described below in the section "Flux flow resistivity".
In Fig. 4 we report 1/Q 0 (main panel) and the normalized f 0 (inset) as a function of magnetic field, at the fixed temperature of T = 6 K . Q degradation starts at quite low applied fields, while above about 0.5 T, a linear trend for 1/Q 0 (H dc ) sets up. In order to understand whether this degradation can be ascribed to vortex motion, Figure 3. Upper critical field as a function of the reduced temperature. Solid symbols were deduced from CPWR measurements of the critical temperature for different applied DC magnetic fields, with T c0 the critical temperature at zero field (see text). Red crosses were obtained from DC resistance measurements, by the 90% of R n criterion. Open symbols are evaluations of B c2 at T = 6 K , obtained by the TDGL analysis of ρ ff (T) data, as described below . Solid lines are fit of the CPWRs data (see text), extrapolated to lower temperatures (dotted lines) to check the matching with the points from TDGL analysis. www.nature.com/scientificreports/ the field of first vortex penetration should be evaluated. Actually, while L (0) is expected to increase with the Ti amount, the coherence length decreases: according to Ref. 18 , a reasonable estimation for our composition should be ξ ≈ 4.5 nm . This gives a lower critical field at T = 6 K of µ 0 H c1 (6 K) = � 0 [ln(κ) + 1/2]/ 4π 2 L ≈ 6.9 mT , where 0 is the flux quantum and κ = /ξ is the Ginzburg-Landau parameter. The applied perpendicular field at which the first vortices are expected to penetrate the CPWR is 19 is the thickness and w the width of the stripline) which for our CPWRs at T = 6 K is about 0.4 mT. This value is higher than the amplitude of the microwave field, but it turns out that all the resonances whose parameters are reported in Fig. 4, with the exception of the first three points, were measured for H dc > H p , i.e. in the presence of vortices.
Once we established that a source of dissipation and Q degradation for our CPWRs can be ascribed to the presence of vortices, and more precisely to vortex dynamics, we should also consider the existence of vortexpinning forces that act to hamper vortex motion. The main pinning parameters can be deduced by our measurements, as follows.
Vortex-motion complex resistivity and pinning parameters. The microwave behavior of a superconductor can be described through a complex resistivity ρ = 1/(σ 1 − iσ 2 ) , where σ 1 is the quasiparticle conductivity and σ 2 = 1/(µ 0 ω 2 L ) . When vortices penetrate the sample in the mixed state, an additional complex term accounting for vortex motion, ρ vm = ρ vm,1 + iρ vm,2 , should be added, and the effective complex resistivity reads where C is the Campbell penetration depth. The connection to the measured quantities is obtained through the surface impedance Z s , that in the local limit is defined as 20 In the CPWR geometry, the surface resistance and reactance can be expresses as 8  where zero-field measurements yield L and σ 1 , according to the procedure described above. Results for C and σ C in the case of NbTi CPWR at T = 6 K are shown in Fig. 5. The field dependence of C is sublinear, as reported in literature for conventional superconductors 21 .
In order to extract the pinning parameters, the Gittleman-Rosenblum (GR) model is considered, connecting the vortex motion resistivity, as a function of frequency ν , to the viscous drag coefficient η of vortices, their depinning frequency ν p , and the flux-flow resistivity ρ ff 22 : Note that the GR model is in principle valid only in the absence of flux creep. Thus, a proper analysis is needed to exclude significant creep. This was done by means of the dual-frequency DR measurement technique 23,24 . The measurements and corresponding analysis, briefly summarized in the "Methods" section and fully described in a forthcoming paper, allowed us to determine that a relevant role of flux creep can be ruled out in the measured NbTi samples.
Then, before using the GR model-that was initially formulated for nearly straight perpendicular vortices in thin films-one should also check whether the whole vortices or at least significant segments of them are driven by the rf field, otherwise the surface impedance could mainly be controlled by the bending rigidity of the vortex (5) R s = ωL l w eff /Q 0 ; X s = ωL k l w eff = µ 0 ω , www.nature.com/scientificreports/ lattice. In the presently explored regime H/H c1 ≫ 1 , i.e. the vortex density is relatively high and uniform, so that the actual penetration depth of rf fields can be evaluated as ℜ(˜ ) , with ˜ as defined by Coffey and Clem 25 . This length must be compared to the half-thickness of the film, since in CPWRs the rf field penetration occurs from both the surfaces of the film. We obtained ℜ(˜ ) = 0.51 µm and 0.97 µm at the fields of 1 T and 4 T, respectively, for CPWR#1 ( 0.5 µm and 0.83 µm for CPWR#2), to be compared with the half-thickness of 0.7 µm ( 1.2 µm ). Thus, the entire vortex lines or the most part of them are effectively driven by rf currents. Moreover, we made measurements in a regime of low rf current densities, i.e. about 3000 A/cm 2 at 6 K and 1 T. This means that typical displacements during vortex oscillations are 0.1 nm . In such regime, all the pinning centers, both strong and weak, are active along the whole vortex line, and bending effects are reduced. We thus conclude that the GR model is fully applicable in this case. Starting from Eqs. (3) and (8), the parameters ν p and η can be expressed as a function of C and σ C , as Finally, according to the GR model, the pinning constant k p can be calculated as The k p constant, characterizing the pinning force strength, and the vortex viscosity, depending on the relaxation rate of quasiparticles inside the vortex core, are shown in Fig. 6. The reported values indicate rather strong pinning, suggesting that besides grain boundaries, pinning could benefit from the possible presence of smaller and densely distributed defects (e.g. Ti precipitates). In fact, it is known that in samples with Ti concentrations exceeding 45 wt% (61 at.%), as for our CPWR samples (69 at% Ti), the main flux pinning source is from normal conducting α-Ti precipitates 26,27 . Incidentally, the presence of Ti precipitates could also explain other peculiarities, i.e. the measured values of normal state resistivity and upper critical field smaller than values reported in literature for the nominal composition, since the presence of Ti precipitates would imply a lower Ti content in the alloy itself, and thus an actual composition different from that overall measured by EDX on a macroscopic area. Figure 6 also shows virtually absent (a) or moderate (b) sample-to-sample variability.
In the comparison with the DR data, one should take into account that DR measurements were done on NbTi films with a different nominal composition (62.3 at% Ti, i.e. lower than CPWR samples and at the limit for the formation of Ti precipitates) and therefore they could present differences in their pinning properties, as well as they actually show different normal-state resistivity values. In Fig. 7a we report the depinning frequency ν p for the two CPWRs (again showing reasonable sample-to-sample variability), compared to the results of the DR measurements. DR data show lower ν p , reflecting a lower pinning force, even by including the different values of the normal state resistivities, and this is consistent with the hypothesis that precipitation of normal conducting α-Ti could be negligible for this Ti concentration. Nevertheless, the order of magnitude of ν p is the same for CPWR and DR, and also the slope of the field dependence is very similar, suggesting the existence of a collective pinning regime in all the samples. Where CPWR and DR data collapse is in the scaling of the flux flow resistivity, shown in Fig. 7b. In fact, despite a significant difference in the normal state properties, the curves of the flux Flux flow resistivity. Figure 8 reproposes the flux-flow resistivity normalized to the normal-state resistivity ρ ff /ρ n as a function of the applied magnetic field, for the data collected at the fixed temperature of T = 6 K . The ρ n (T c ) = 22.5 µ� cm value was determined from the four-contact DC resistivity measurements of an Hall bar from the same NbTi films, as described in the "Methods" section. As a first attempt, this plot should be discussed in the framework of the Bardeen-Stephen (BS) model 28 . According to BS, in dirty s-wave superconductors the dissipation due to vortex motion is to be ascribed to the quasiparticle currents flowing in the normal cores of vortices, resulting in ρ ff /ρ n = B/B c2 . Clearly our data do not follow this trend, since they show an upward curvature. The reasons for this behavior must be found beyond the elementary theories of flux flow. Deviations  www.nature.com/scientificreports/ from the simple proportionality ρ ff ∝ B (i.e. ρ ff simply proportional to the density of vortices) emerge when the suppression of the order parameter with increasing field and modifications in the current patterns around vortices are considered. Among the most convincing approaches, the time dependent Ginzburg-Landau (TDGL) theory can be used to determine the detailed ρ ff (T, B) behavior beyond the BS approximation. According to the applicable mean-field result 29 , To the first order α ≈ 0.4 , remarkably independent of temperature, field, and material parameters. The field dependence is about linear only when B ≪ µ 0 H c2 , where ρ ff /ρ n ≃ αB/B c2 , but even in this case the BS result is not recovered since α < 1 . At increasing B an upward curvature is expected.
For each CPWR we found the B c2 (6 K) value that adapts the Eq. (12) curve to the experimental data, resulting in α values close to the expected 0.4. A rather good agreement was found for α = 0.35−0.36 for the two CPWRs. Moreover, the B c2 (6 K) values obtained by this method were plotted above in Fig. 3, to show their consistency with the B c2 values evaluated from direct measurements at higher temperatures.
Comparing CPWR with DR-measured samples, the overlapping ρ ff /ρ n reported in Fig. 7b implies that the ratios α/B c2 in all the samples (in the different field range explored) are essentially the same. Assuming the expected universal value α = 0.4 also for the DR sample, this would imply compatible values of B c2 between the CPWR and DR samples also considering the uncertainty in the trend extrapolation to higher fields for the DR measurements, more limited in terms of the maximum H dc ≤ 1.2 T applied.
Overall, the flux flow resistivity analysis shows that data can be understood within the frame of the TDGL theory. In fact, a TDGL model approach was recently used to explain the vortex-induced nonlinear rf response of Nb SRF cavities 30,31 .

Conclusions and perspectives
The adopted coplanar-waveguide-resonator technique revealed to be reliable to study the high-frequency properties of superconducting films in the presence of DC magnetic fields. By analyzing a complete set of measurements, we were able to obtain the absolute value of the penetration depth, the complex impedance and the contribution of the vortex motion to the complex resistivity, that was one of the main targets of the work. Moreover, to access the pinning parameters we analyzed the complex impedance within the formalism of the Campbell penetration depth. The complete parameter set shown in Figs. 6 and 7 comprehensively characterizes the vortex-pinning properties of the material, with moderate sample-to-sample variability.
Measurements in the GHz frequency range offer the advantage to investigate both the pinning properties and the flux flow resistivity. The last quantity was discussed and general scaling laws with normalized field were investigated. It turns out that the normalized flux flow resistivity follows the trend predicted by the time dependent Ginzburg-Landau theory remarkably well. This final result validates the whole CPWR measurement and analysis process, thus promoting the CPWR-based approach as a promising technique for the study of superconducting films in view of applications in high frequency and high field conditions.
One of the merits of the adopted approach is that all the quantities needed for the analysis, included the normalization factors, were extracted from the measurements of the very same films and devices, thus avoiding the use of data from literature which is problematic due to the high variability of parameters with the Ti content and preparation process of the samples. In fact, we could compare the CPWR results to a different microwave analysis (based on the use of a dielectric-loaded resonator), but only for NbTi films with a different composition, even if deposited in the same sputtering apparatus. Notwithstanding differences in the pinning force, the depinning frequency shows a quite reasonable agreement, and the normalized flux flow resitivities, as deduced by the two techniques, remarkably overlap.
In conclusion, the ability to arrange the whole data-set in a consistent framework of theory and pinning models is the real added value of this work, since it defines an analysis procedure and it constitutes the base to predict the behavior of the material under the application conditions and thus to design practical devices.
In perspective, the CPWR technique offers interesting possibilities for the investigation of nonlinear effects due to vortex motion, since in this configuration (rf current peaks at the edges) the nonlinearity threshold is easily achievable even with the moderate rf power of vector network analyzers. We are currently performing these studies on NbTi films and other Nb-based compounds.

Methods
Thin film deposition. The deposition of the NbTi films onto quartz substrates was carried out by DC magnetron sputtering using a 4″ commercial planar source with a target-sample distance of 11 cm. We deposited two sets of films, with different experimental conditions: the chamber was baked at 600 °C for 24 h and 48 h for the two runs, giving a base pressure before the process at 550 °C of 1 · 10 −7 mbar and 2 · 10 −8 mbar , respectively. The processes were then carried out at 550 • C fixing the current at 1.5 A (corresponding to a current density of 0.148 A/cm 2 and to a voltage of about 400 V). Argon was used as working gas, with a pressure of 6 · 10 −3 mbar . The deposition times were 20 min in the first process and 40 min in the second one. After deposition and cooling to room temperature, the pressures of the system were 2 · 10 −8 mbar and 7 · 10 −10 mbar , respectively. The thicknesses of the films were measured using a micro-profilometer and checked by Atomic Force Microscopy (AFM) on the devices, and resulted to be 1.4 μm for the first run samples and 2.4 μm for the second run (more similar to axion cavities coating). www.nature.com/scientificreports/ Device fabrication. Coplanar waveguide resonators were produced through a combination of classical lithography, sputtering depositions and reactive ion etching. Two devices were analyzed in this work (labelled as #1 and #2), representative of the two deposition processes described in the previous sub-section. The patterning was obtained by means of an optical lithographic laser system with a resolution of 0.9 μm. The photoresist used during the design patterning did not provide an adequate selectivity with the NbTi film, therefore the polymer mask was not durable enough to guarantee a complete etching of the superconductive film. In order to overcome this issue, an approach based on a hard mask was adopted, in particular depositing 120 nm of aluminium sputtered in UHV conditions ( 2 · 10 −8 mbar) through an argon plasma ( 10 −3 mbar) with a deposition rate of 1.3 nm/s. The Al layer in excess was removed in an acetone bath, useful for the photoresist dissolution.
In these conditions, a highly directional dry etching was performed, combining the chemical and physical resistance of the aluminum to the fluorine-based gas mixture for incompatibility, and to the plasma etching exploiting its thickness, respectively. Reactive Ion Etching relied on SF 6 at 2 · 10 −2 mbar, glowed in plasma through an ICP power of 800 W and RF power of 50 W, leading to an etching rate of 2 nm/s. The hard mask was then removed using a wet etching in a NaOH bath.
Ancillary measurements. The morphology of the films was analyzed by field emission scanning electron microscopy (FESEM), showing the presence of grains with a mean size of about 200 nm. FESEM analysis was finally performed after the device fabrication, to check the high quality of the transferred pattern to the NbTi film, with particular attention to the sharpness of the edges.
The composition of the films was measured by EDX (Energy Dispersive X-ray Spectroscopy) over an area much greater than the dimension of the grains, and turned out to be Nb 0.31 Ti 0.69 (atomic percentage). This composition was expected to maximize H c2 , rather then T c 32 .
DC resistivity measurements were performed to find out parameters useful for the microwave analysis, such as the low-temperature normal state resistivity. Hall bars 15 µm in width and 1 mm in length were produced with the same method described above and their resistivity was measured by a standard four probe current-biased technique. To cancel the thermoelectric voltage offset and its possible time drift, the delta method was applied: each data point is the average of three voltage readings carried out alternating the polarity of the bias current 33 . In addition, to avoid sample heating, each current pulse was 200 ms long with an interval of 5 s between two measurements. It came out that the normal state resistivity at T T c is ρ n = 22.5 µ� cm , while a residualresistivity ratio of 1.74 was found.
CPWR microwave measurements. The CPWR was capacitively coupled to the readout circuit in a brass package in tight thermal contact with the cold finger of a closed cycle cryocooler. The complex transmission coefficient S 21 (ratio of the voltage transmitted to the incident voltage), as a function of the driving frequency, f, was detected by a vector network analyzer in different conditions of temperature, DC magnetic field (applied perpendicular to the film plane), and rf input power. The resonance frequency f 0 and the loaded quality factor Q L were obtained from a fit of the experimental magnitude of the transmission coefficient: and the unloaded quality factor, Q 0 , was then obtained accounting for the coupling coefficients. Figure 1 shows examples of such resonance curves measured at different temperatures for the NbTi CPWR #1. The temperature dependence of the resonance frequency and of the inverse of the unloaded quality factor (obtained starting from the fit of the single curves by Eq. 13), is shown in the inset. These parameters were used to determine the London penetration depth and the surface resistance, as detailed above.
Dielectric-loaded resonator technique. Another NbTi sample was prepared to be studied, without patterning, through surface impedance measurements by means of a dielectric-loaded cylindrical resonator (DR) 6 . The sample used with the DR is ≈ 1.7 μm thick (evaluation from calibrated deposition rate and deposition time), with 62.3 at.% Ti (46% weight), grown on a square quartz substrate, 1.2 mm thick and with an area of 25 × 25 mm 2 . The normal-state resistivity at T T c , as obtained by the van der Pauw method with the nominal thickness, is 54 µ cm.
Following a perturbation approach, the sample was placed on one of the resonator bases and covered with a thin metal mask with a circular hole (diameter = 17mm) to preserve the circular symmetry. The sample thus contributes to the resonator unloaded quality factor Q 0 and resonance frequency f 0 . From the Q 0 and f 0 field induced variations, the sample surface impedance variations with H dc can be obtained 6 : The resonator can be excited at two electromagnetic modes (transverse electric modes TE 011 and TE 021 with resonant frequencies f 01 = 16.5 GHz and f 02 = 26.7 GHz, respectively), so that �Z s (H) (and the extracted ρ vm ) is measured at the two distinct frequencies f 01 and f 02 24 . The availability of four observables (real and imaginary parts of the surface impedance at two distinct frequencies) makes possible to use more extended ρ vm models 23,25,34 with respect to the GR model and to check its basic assumption about the negligibility of creep. Measurements www.nature.com/scientificreports/ were performed at fixed, selected temperatures between 5 K and T c ≃ 9.1 K, by applying H dc ≤ 1.2 T normally to the sample surface. Similarly to what observed with the CPWR, the penetration depth of the microwave e.m. fields in the NbTi film varies, changing the regime from bulk to thin film, by increasing T and H dc . The finite thickness effects are properly taken into account resorting to surface impedance transformation relations 6,20 . It can be observed that at 6 K these effects impact the vortex parameters values, with respect to the bulk limit, by ∼ 10% . Vortex parameters can be extracted resorting to the Coffey-Clem (CC) model 25 : where χ ∈ [0, 1] as an adimensional thermal creep factor and ν c is a characteristic frequency, which reverts to ν p for χ → 0 . Indeed, the CC model expression for ρ vm with zero creep is equal to the GR expression, Eq. (8). The full analysis of the surface impedance data (to be published in a forthcoming paper) allowed us to assess that creep effects at low T (6 K included) are indeed negligible, confirming the approach taken in analysing CPWR measurements. The other two vortex parameters, ρ ff and ν p , are shown in Fig. 7.